Evanescent channels and scattering in cylindrical nanowire heterostructures 



P. N. RacecQ 

Weierstrafi-Institut filr Angewandte Analysis und Stochastik, Mohrenstr. 39 10117 Berlin, Germany and 
National Institute of Materials Physics, PO Box MG-7, 077125 Bucharest Magurele, Romania 

E. R. RacecQ 

Institut filr Physik, Technische Universitdt Cottbus, 
Postfach 101344, 03013 Cottbus, Germany and 
Faculty of Physics, University of Bucharest, PO Box MG-11, 077125 Bucharest Magurele, Romania 

H. NeidhardlE] 

Weierstrafi-Institut fiir Angewandte Analysis und Stochastik, Mohrenstr. 39 10117 Berlin, Germany 

We investigate the scattering phenomena produced by a general finite-range nonseparable poten- 
tial in a multi-channel two-probe cylindrical nanowire heterostructure. The multi-channel current 
scattering matrix is efficiently computed using the R-matrix formalism extended for cylindrical co- 
ordinates. Considering the contribution of the evanescent channels to the scattering matrix, we 
are able to put in evidence the specific dips in the tunneling coefficient in the case of an attractive 
potential. The cylindrical symmetry cancels the "selection rules" known for Cartesian coordinates. 
If the attractive potential is superposed over a non-uniform potential along the nanowire, then 
resonant transmission peaks appear. We can characterize them quantitatively through the poles 
of the current scattering matrix. Detailed maps of the localization probability density sustain the 
physical interpretation of the resonances (dips and peaks). Our formalism is applied to a variety of 
model systems like a quantum dot, a core/shell quantum ring or a double-barrier, embedded into 
the nanocylinder. 
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I. INTRODUCTION 

In the last years, there is an increased interest in study- 
ing nanowire-based devices due to their broad applica- 
tion area. They can be used as field-effect transistors 
(FET)^ or gate-all-around (GAA) FET^^"^, nanowire 
resonant tunneling diodes^ ^, solar cells as integrated 
power sources for nanoelectronic system^, laser^, and 
also as qubits^. Their structural complexity has progre- 
sively increased, and the material composition includes 
III-V materials but also, so attractive for semiconductor 
industry, group IV materials. 

Description of electrical transport and charge distri- 
bution in nanowire-based devices has to be done quan- 
tum mechanically, and the most appropriate method for 
such open systems is the scattering theory. Due to the 
confinement of the motion inside the nanowire, the elec- 
trons are free to move only along the nanowire direc- 
tion, so that these systems are also called quasi-one- 
dimensional systems. Since many nanowires have circular 
cross-sectional shape^ ^ ^, we present in this work a gen- 
eral method, valid within the effective mass approxima- 
tion, for solving the three-dimensional (3D) Schrodinger 
equation with scattering boundary conditions in cylindri- 
cal geometries. The azimuthal symmetry suggests to use 
cylindrical coordinates, with z axis along the nanowire. 
This reduces the scattering problem to two dimensions: 
r and z directions. Its solut ion is found numeric ally 
using the R-matrix iorm^W^m^^^^^^^^^^ ex- 



tended for cylindrical coordinates. 

An interesting effect in a multi-channel scattering 
problem is that as soon as the potential is not any- 
more separable, the channels get mixed. If further- 
more the scattering potential is attractive, then it leads 
to unusual scattering properties, like resonant dips in 
the transmission coefficient just below the next chan- 
nel minimum energy. As it was shown analytically for 
a 5 scattering potentiali^ or for point scatterer^, and 
later on for a finite-range scattering potential^^^ the 
dips are due to quasi-bound-states splitting off from a 
higher evanescent channel. So that evanescent chan- 
nels can not be neglected when analyzing scattering in 
two- or three-dimensional quantum systems. These find- 
ings were recently confirmed numerically for a Gaussian- 
type scatterei^^ and also for a quantum dot or a quan- 
tum ring^ embedded inside nanowires tailored in two- 
dimensional electron gas (2DEG). 

It is the aim of this work to show that we could find 
the same features in the case of a cylindrical nanowire. 
Furthermore, in cylindrical nanowires, due to the three- 
dimensional (3D) modelling, every magnetic quantum 
number defines a two-dimensional (2D) scattering prob- 
lem, with different structure of dips for the same scatter- 
ing potential. Also, the cylindrical symmetry does not 
forbid any intersubband transmission, so that we could 
find dips in front of every plateau in the transmission co- 
efficient. We apply our method to a variety of model sys- 
tems like quantum dot, quantum ring or double-barrier 
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heterostructure, embedded inside the nanowire. 



II. MODEL 

We consider a cylindrical nanowire with a constant po- 
tential on the surface. Inside the wire the electrons are 
scattered by a potential of finite extend. 



A. Scattering problem for the cylindrical geometry 

In the effective mass approximation, the envelope func- 
tion associated to the energy E satisfies a Schrodinger- 
type equation 



V{r) 



«'(r) = ^*(r). 



(1) 



We use the symbol /i to denote the effective mass of 
the electrons, while m will denote the magnetic quan- 
tum number. As long as there are not split gates on 
the surface of the nanowire, the potential energy V{r) is 
rotational invariant 



V(r) = V{r,z) 



(2) 



and nonseparable in a small region of the structure called 
scattering region. The z-axis was considered along the 
nanowire as shown in Fig. [l] 
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FIG. 1: The geometry of the 2D scattering problem. 

A scattering potential which does not explicitly depend 
on the azimuthal angle imposes the eigenfunctions of 
the orbital angular momentum operator as solutions 
ofEq. ^ 



(3) 



where m = 0,±1,±2,... is the magnetic quantum num- 
ber. This is an integer number due to the require- 



ment that the function 



should be single- valued. 



The allowed values of the energy E and the functions 
i^m^E] r, z) are determined from the equation 



2/i \dr'^ r dr 



Eil;m{E;r, z). 



(4) 



Due to the localized character of the scattering potential 
it is appropriate to solve Eq. Q within the scattering 
theory. In such a way, every magnetic quantum num- 
ber m defines a two-dimensional (2D) scattering prob- 
lem. Furthermore, these 2D scattering problems can be 
solved separately if the scattering potential is rotational 
invariant. How many of these problems have to be solved, 
depends on the specific physical quantity which has to be 
computed. 

The potential energy which appears in Eq. Q has 
generally two components: 



V{r,z) = V^{r) + Vseatt{r,z). 



(5) 



The first one, V±(r), describes the lateral confinement of 
the electrons inside a cylinder of radius R and is trans- 
lation invariant along the nanowire. We consider a hard 
wall potential 



V^{r) 



0, 0<r < R 
oo, r > R 



(6) 



suitable for modelling either free-standing nanowires or 
nanowire transistors with no gate leakage current. Both 
situations correspond to the state-of-the-art devices. 

The scattering potential energy inside the nanowire, 
Vscatt{r, z), has generally a nonseparable character in a 
domain of finite-range and is constant outside this do- 
main. We consider here the nonseparable potential lo- 
calized within the volume defined by the boundaries ±dz 
and i?, see Fig. [l] 



Vscatt{r,z) = < 



Vi, re%R],z<-d, 

Vs{r,z), re [0,R],-d,<z<d, . (7) 

^2, re [0,R],z>d, 



There are not material definitions for the planes z = :^dz . 
Usually, they are chosen inside the highly doped regions 
of the nanowire characterized by a slowly z-varying po- 
tential, practically by a constant potential. These regions 
play the role of the source and drain contacts. 



B. Scattering states 

In the asymptotic regions, \z\ > dz i.e. source 
and drain contacts, the potential energy is separable in 
the confinement and the transport direction (F(r, z) = 
V±{r)^Vs) andEq. Q can be directly solved using the 
separation of variables method 

^m{E;r,z)=^{rMz). (8) 

The function (j){r) satisfies the radial equation 



d^ 



1 d 



2/i dr^ r dr 



m 



{r) = E^<P{r), (9) 
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while (f{z) satisfies the one-dimensional Schrodinger type 
equation 



^{z) = {E - E^Mz), (10) 

and 



where 5 = 1 stays for the source contact {z < —dz^ 
5 = 2 for the drain contact {z > dz). 

Due to the electron confinement inside the cylinder of 
radius the solutions of Eq. ([9| are given in terms of 
the Bessel functions of the first kind, J^, 



)(r) = 



V2 



71 



(11) 

where Xmn is the nth root of Jm(^)- The eigenfunctions 
(t)^\r)^ called transversal modes^ form an orthonormal 
and complete system of functions. The corresponding 
eigenenergies are 



E 



(m) _ 



2/1 \ R 



1,2, 



(12) 



and they depend only on the effective mass and the radius 
of the cylindrical nanowire. It is worth to mention here 



that l^n and and £^Y^ depend only on |m|. 

Every transversal mode together with the associated 
motion on the transport direction defines the scattering 
channel on each side of the scattering area (5 = 1 for 
the source contact and 5 = 2 for the drain contact). The 
scattering channels are indexed by {mns) for each E. 

If the total energy E and the lateral eigenenergy E^^ 
are fixed, i.e. for every E, m and n, there are at most 
two linearly independent solutions of Eq. (10). In the 



asymptotic region they are given as a linear combination 
of exponential functions 



r^sn 



(^) 



iki 



-iki 



Z < 



(13) 



where Ag, Bg, Cg and Dg are complex coefficients de- 
pending on n, m and E for each value of s = 1,2. The 
wave vector is defined for each scattering channel {mns) 
as 



{E-E 



(m) 



(14) 



where ko = 7r/2dz and uq - 
conducting or open channels 



h^kQ/2ii. In the case of 



E-E 



(m) 



>o. 



(15) 



ksnm are positive real numbers and correspond to propa- 
gating plane waves. For the evanescent or closed channels 



E-E 



(m) 



<o. 



(16) 



ksnm are given from the first branch of the complex 
square root function, ksnm = ^l^snml^ and describe ex- 
ponentialy decaying functions away from the scattering 
region. Thus, the number of the conducting channels, 
Nsni{E)^ 5 = 1, 2, m > 0, is a function of energy, and 
for a fixed energy E this is the largest value of n, which 
satisfies the inequality (15) for given values of s and m. 



Each conducting channel corresponds to one degree of 
freedom for the electron motion through the nanowire 
and, consequently, there exists only one independent so- 
lution of Eq. (|4| for a fixed channel {mns) associated 
with the energy E^ ipnm{E]r^ z). For describing further 
on the transport phenomena in the frame of the scatter- 
ing theory it is convenient to consider this solution as 
a scattering state^ i.e. as a sum of an incoming compo- 
nent on the channel {mns) and a linear combination of 
outgoing components on each scattering channel. In a 
convenient form^^, the scattering function incident from 
the source contact (s = 1) is written as 



V'ill(i?;r,z) 



e{Ni^{E) - n) 



27r 



and the scattering function incident from the drain contact {s = 2) as 



^^nl{E;r,z) 



0{N2m{E) - n) 



27r 



> ^2 



(17a) 



n' = l 

e 



(17b) 



The step functions in the above expressions, with 0{x > 0) = 1 and 0{x < 0) = 0, assure that the scattering 
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functions are defined only for conducting channels. 

The three-dimensional scattering states, solutions of 
Eq. ([T]) for rotational invariant geometries can be now 
written as 



Being eigenfunctions of an open system, they are or- 
thonormalized in the general sense^ 



/27r 



(18) 



/oo pR />27r 

dz drr de<i>mE;r,e,z)^l:^,l{E';r,e,zy 
-OO ^0 ^0 



/ 6 oaf Sfi 



-'mm' '-'ss' ^nn 



5{E - E') 

Hsnm 



r 



(19) 



where gsnm{E) is the ID density of states, gsnm{E) = 

IJi/[h^ksnm{E)]. 

The physical interpretation of the expressions ( pT] ) is 
that, due to the nonseparable character of the scattering 
potential, a plane wave incident onto the scattering do- 
main is reflected on every channel - open or closed for 
transport - on the same side of the system and trans- 
mitted on every channel - open or closed for transport - 
on the other side. The reflection and transmission am- 
plitudes are described by the complex coefficients S^^) 

and S^^, with s ^ s\ respectively, and all of them 
should be nonzero. These coefficients define a matrix 
with Nim{E) -\- N2m{E) infinite columns. For an elegant 
solution of the scattering problem we extend S^'^^E) 
to an infinite square matrix and set at zero the ma- 
trix elements without physical meaning, S^T^, = 0, 
n > Nsm{E), 5 = 1,2. In this way we define the wave 
transmission matrix^^ or the generalized scattering ma- 
trioS^. This is not the well-known scattering matrix (cur- 
rent transmission matrix) whose unitarity reflects the 
current conservation. The generalized scattering matrix 
is a non-unitary matrix, which has the big advantage that 
it allows for a description of the scattering processes not 
only in the asymptotic region but also inside the scatter- 
ing area. 

For the sake of simplicity and also considering that for 
rotational invariant potentials the 2D scattering prob- 
lems generated by every magnetic quantum number m 
can be solved separately, the index m will be omitted in 
the following subsections. 



C. R-matrix formalism for cylindrical geometry 

The scattering functions inside the scattering region 
are determinated using the R-matrix formalism, i.e. 
they are expressed in terms of the eigenfunctions cor- 
respond ing to the closed c ounterpart of the scattering 
prQ]QlejjPLll i | i2 | i3 | i4 | i5 | i 6lT7l^ opinion this is a more 

appropriate method than the common mode space ap- 
proach which implies the expansion of the scattering 
functions inside the scattering area in the basis of the 



transversal modes ^n(^)- As it is shown in Ref. |22|27j 
the mode space approach has limitations for structures 
with abrupt changes in the potential or sudden spatial 
variations in the widths of the wire; it breaks even down 
for coupling operators that are not scalar potentials, like 
in the case of an external magnetic field. In the R- 
matrix formalism the used basis contains all the infor- 
mation about the scattering potential, and this type of 
difficulties can not appear. 

Thus, the scattering functions inside the scattering re- 
gion are given as 



^l:\E;r,z) = Y.a\:j{E)xi{r,z), 



(20) 



with r G [0, i?] and z G [—dz^dz]. 

The so-called Wigner-Eisenb ud fu nctions, xi{^^^)^ 
firstly used in the nuclear physic^^^^, satisfy the same 
equation as ipn\r^ z)^ Eq. (4|, but with different bound- 
ary conditions in the transport direction: Since the 
scattering function ifju^ (r, z) satisfies energy dependent 
boundary conditions derived from Eq. (17) due to 



the continuity of the scattering function and its deriva- 
tive at z = -^dz , the Wigner-Eisenbud function xi (^^ ^) 
has to satisfy Neumann boundary conditions at the 
interfaces between the scattering region and contacts, 
dxi/ dz\^^j^^ = 0, / > 1. The infinite potential outside 
the nanowire requires Dirichlet boundary condition on 
the cylinder surface for the both functions, ipn^ (i?, 2:) = 
and xi{R->^) = 0- The functions xu ^ ^ built a basis 
which verifies the orthogonality- and the closure- rela- 
tion. The corresponding eigenenergies to xi ^^'^ denoted 
by El and are called Wigner-Eisenbud energies. Since the 
Wigner-Eisenbud problem is defined on a closed volume 
with self-adjoint boundary conditions, the eigenfunctions 
Xi and the eigenenergies Ei can be choosen as real quan- 
tities. The Wigner-Eisenbud problem is, thus, the closed 
counterpart of the scattering problem. 

In the case of the one-dimensional system it was re- 
cently proven mathematically rigorous that the R-matrix 
formalism allows for a proper expansion of the scatter- 
ing matrix on the real energy axi^^. In this section we 
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present an extension of the R-matrix formalism for 2D 
scattering problem with cylindrical symmetry. 

To calculate the expansion coefficients (E) we mul- 
tiply Eq. Q by Xi{^^ ^) the equation satisfied by the 
Wigner-Eisenbud functions by ilJn\E;r^ z). The differ- 
ence between the resulting equations is integrated over 
[— (i^, dz] X [0,i?], and one obtains on the right hand 



side the coefficient a\^{E). After an integration by 
parts in the kinetic energy term and using the bound- 



ary conditions one finds and feed in it into Eq. (|20 ) . 



So, the scattering functions inside the scattering region 
{z e [—dz^dz]^ r e [O^R]) are obtained in terms of their 
derivatives at the edges of this domain, 



^ Jo 



dr' r' 



R{E,r\ -dz,r,z) 



d^^:\Ey,z') 



dz' 



R{E 



, _ ^ d^i'\E,r\z^) 
r ,dz,r,z) - 



dz' 



(21) 



where the R-function is defined as 

Xi{r,z)xi{r' ,z') tt 



R{E,r,z,r',z') 



2^ ^ E-Ei 2d, 



(22) 



The functions dipn^ /dz ai z = :^dz are calculated from 
the asymptotic form ( [l7| ) based on the continuity condi- 
tions for the derivatives of the scattering functions on the 
interfaces between the scattering region and contacts. 

With these results the scattering functions inside the 
scattering domain are expressed in terms of the wave 
transmission matrix S 



"^{E'.r.z) 



e{E)[l-S^{E)K{E)R{E;r,z)], (23) 



where the component {sn) of the vector ^ is the scatter- 
ing function ipn\E; r, z),n>l, 5 = 1,2 and S-^ denotes 
the matrix transpose. The diagonal matrix K has on its 
diagonal the wave vectors (p!4|) of each scattering channel 



n, > 1, 5, 5^ = 1, 2, and the vector R{E] r, z) as 

^0 ^Xi{r,z)xi 



R{E',r,z) = 



E-Ei 



where xi is a vector with the components 
1 

{Xi)sn = / xi{rd-'^ydz)(l)n{r)rdr, 

Vko Jo 



(24) 



(25) 



(26) 



n > l,s = 1,2. The diagonal 0-matrix, Ssn,s'n'{E) = 
6{Ns{E) — n) 5ss' ^nn'i n > 1, s = 1,2, assures non-zero 
values only for the scattering functions corresponding to 
the conducting channels. 

Using further the continuity of the scattering func- 
tions on the surface of the scattering area and expanding 



R{E] :^dz^ r) in the basis {(/>n(^)}n>i we find the relation 
between matrixes S and R 

S{E) = [l - 2 (1 + iIl{E)K{E)y^] S{E), (27) 



with the R-matrix given by means of a dyadic product 



1=1 



E-E, 



(28) 



According to the above relation, R is an infinite- 
dimensional symmetrical real matrix. The above form 
allows for a very efficient numerical implementation for 
computing the R-matrix. 

The expression (27) of the S-matrix in terms of the R- 



matrix is the key relation for solving 2D scattering prob- 
lems using only the eigenfunctions and the eigenenergies 
of the closed system. On the base of Eq. (27) the wave 



transmission matrix is calculated and after that the scat- 
tering functions in each point of the system are obtained 
using Eqs. (17) and (23). 



To come back to the dependence on m, we point out 
that the Wigner-Eisenbud functions and energies are m- 
dependent, so that the matrixes R, K and and the 
vector ^(£^;r, z) are m-dependent in relations (27) and 
([23l). 



D. Reflection and transmission coefficients 



Using the density current operator 

j(r) = A (M,(r) VM/(r)* - M/(r)* V*(r)) , (29) 



one can define, as usually, the transmission and reflection 
probabilities.^^ ^(r)* denotes the complex conjugate of 
the scattering wave function (18). 



The r-component of the density current jr{r^0^z) is 
zero in leads, because (pni^) are real functions. The com- 
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ponent of the incident density current is m dependent, with the infinite dimensional matrix 



{jinc{r,e,z))e 



1 1 



/i (27r)2 r 



m|0n(r)|^ 



so that if one sums over ah m values, then they cancel 
each other. This is also valid for the reflected and trans- 
mitted current fluxes. What remains is the 2:-component 
of the particle density current jzi^, 0, z), which integrated 
over the cross section of the nanocylinder with the cor- 
responding measure, rdr, provides the very well known 
relations for the transmission and reflection probabili- 
ties. The probability for an electron incident from source, 
5 = 1, on channel n to be reflected back into source on 
channel n' is 



^ ^In^int |2 
-^nn' I, \^ln,ln'\ •> 



(30) 



and the probability to be tranmitted into drain, 5 = 2, 
on channel n' is 



7^(1) 

nn' 



^2ri/\Qt 1 2 

7 l^ln,2nM • 



(31) 



The reflection and transmission probabilities for evanes- 
cent (closed) channels are zero. The total transmission 
and reflection coefficients for an electron incident from 
reservoir 5 = 1 are defined as 



(32) 



More detailed properties of the many-channel tunneling 
and reflection probabilities are given in Ref. [30 , but 
note that our indexes are interchanged with respect to 
the definitions used there. Of course, all these coefficients 
are m-dependent. 



E. Current scattering matrix 

Further, we define the energy dependent current scat- 
tering matrix as 

s{E) = K^/2(£;)0(£;)s(£;)K-i/2(£;), (33) 

so that its elements give directly the reflection and trans- 
mission probabilities 



(2) 



(34) 



(35) 

The diagonal 0-matrix assures that the matrix ele- 
ments of S are nonzero only for the conducting channels, 
for which the transmitted flux is nonzero. Using the R- 
matrix representation of S, Eq. (27), we find from the 
above relation 

S{E) = S{E) [1 - 2(1 + in{E))-^] S{E), (36) 



(2)/ 



n{E) =uoJ2 1^ = Ki/2(i?)R(i?)Ki/2(i?) (37) 



1 = 1 



and the column vector 



(38) 



with / > 1. According to the definition (37) the matrix 
is a symmetrical one, ft = 11-^, and from Eq. (36) 
it follows that S also has this property, S = . On 
this basis one can demonstrate that the tunneling coeffi- 
cient characterizes one pair of open channels irrespective 

of the origin of the incident flux T^^^, 

2 



J2n'ln 



Sin. 



2n' 



This is a well-known property of the 

transmission through a scattering system and it shows 
that the current scattering matrix used here is properly 
defined. The restriction of S-matrix to the ope n chan- 
nels is the well known current scattering matrij^* -*^^ * -*^^ *, 
commonly used in the Landauer-Biittiker formalism. For 
a given energy E this is a {Ni -\- N2) x (Ni + N2) matrix 
which has to satisfy the unitarity condition, according to 
the flux conservation. 



The relation ( 36 ) is the starting point for a resonance 
theory-*^^'^^, which allows for an explicit analytical expres- 
sion for the transmission peak as a Fano resonance with 
a complex asymmetry parameter. 

In the numerical computations, the matrixes S, R, 12, 
S and have the dimension 2N x 27V, and the vectors 
Xh have 2N components, where N is the number 

of scattering channels (open and closed) taken numeri- 
cally into account. The number of the Wigner-Eisenbud 
functions and energies computed numerically establishes 
the maximum value for the index 



III. CYLINDRICAL NANOWIRE 
HETEROSTRUCTURE MODEL SYSTEMS 



The formalism presented above is general and can 
be applied to a variety of the cylindrical nanowire het- 
erostructures. We consider a series of heterostructures 
embedded in an infinite cylindrical nanowire of radius 
R = 5nm and effective mass /i = 0.19mo (corresponding 
to transverse mass in Silicon). We set in all our compu- 
tations dz = 16nm, see Fig. [l] and the total number of 
channels (open and closed) N = 8. In our calculations, 
the results do not change if more channels are added. In 
Fig. [2] are plotted the energies of the transversal modes 
^j^n u^^il leF, for different magnetic quantum numbers 



m, according to (12). The difference between two suc- 



cessive energies of the transversal modes is m-dependent, 
due to the roots Xmn of the Bessel functions Jm- 
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FIG. 2: (color online) Energies E^_^l of the transversal modes 
for a cylinder with R — 5nm and fi — 0.19mo. 



A. Quantum dot embedded into the nanocylinder 



1. Same radius as the host cylinder 



In Fig. 3(a) is sketched a cyHndrical quantum dot em- 
bedded into a cyHndrical nanowire with the same radius. 
This kind of structures and even compositionally modu- 
lated, called also "nanowire superlattice"^^, are already 
realized technologically on different materials basis, as is 
summarized in a recent review article^^. 

Depending on the band-offsets between the dot mate- 
rial and the host material the potential produced by the 
dot can be repulsive, yielding a quantum barrier, or at- 
tractive, yielding a quantum well. As it is mentioned in 
Ref. |3Ti the interfaces between the dot and the host ma- 
terial may be considered sharp for nanowires with diam- 
eter less than 20nm. We consider here that the dot yields 
an attractive potential V(r, z), represented in Fig. |3(b)| 
by a rectangular quantum well of depth Wi) = —0.5eV 
and width w = 8nm. 

The total tunneling coefficient T^^^ versus the incident 
energy E is plotted in Fig. |4j for different magnetic quan- 
tum numbers m and different quantum well depths W^. 
In the absence of the quantum well, = O.OeF, one can 
recognize the abrupt steps^^ in the tunneling coefficient. 
The transmission increases with a unity, every time a 
new channel E^_^^ becomes available for transport, i.e. 
becomes open. The length of the plateaus is given by 
the difference between two successive transversal mode 
energies, which differ for different m values. Due to the 
almost square dependence of the transverse energy lev- 
els on the channel number, the length of the plateaus 
increases. Increasing the depth of the well, deviations 
from the step-like transmission appear. There are no ef- 
fects due to the influence of the evanescent channels^^, 
because the scattering potential for this configuration re- 



Source 



] 





Drain 




H 







w 

(a) 




FIG. 3: (color online) (a) Sketch of a cylindrical quantum dot 
embedded into a nanowire with the same radius. We consider 
that the dot yields an attractive potential V{r, z), represented 
in (b) by a rectangular quantum well of depth Wb = —0.5eV 
and width w — 8nm. 



mains furthermore separable in the confinement and the 
transport direction, V{r,z) = U{r) -\- W{z). 

The spectral representation of the 2D Hamiltonian in 
this situation is a superposition of the spectrum of each 
channel, without being perturbed by channel mixing. 
Considering an attractive potential in z-direction, there 
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FIG. 4: (color online) Total tunneling coefficient as a function 
of inci dent e nergy E for the scattering potential represented 
in Fig. 3(b) for different magnetic quantum numbers m and 
various values of the well depth Wh — O.Oey (continuous line), 
Wh = 0.05ey (dashed line), and Wh = CSeV (dot-dashed 
line) . 
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is always at least one bound stat^^^I^ below the c on- 
tinuum spectrum for every channel n, see Fig. |7(a)[ a). 
In turn, the bound states of the higher channels n get 
embedded in the continuous part of the lower channels, 
forming bound states in continuum (BIG). Since the po- 
tential is separable, there is no mix of states, and the 
BIG states can not be seen as scattering states. 



2. Surrounded by host material 

Further we study a cylindrical dot embedded into the 
nanocylinder, but whose radius is smaller than the 
cylinder radius so t hat th e dot gets surrounded by the 

We consider here again the 



host material, see Fig. 5(a) 



case that the dot yields an attractive potential, i.e. a rect- 
angular quantum well, plotted in Fig. |5(b)[ Even we have 
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FIG. 6: (color online) Total tunneling coefficient as a func- 
tion of total incident energy E for different magnetic quantum 
numbers m for a cylindrical dot surrounded by the host mate- 
rial as in Fig. [5] A detailed view around the channel minima 
is presented in the insets. 



FIG. 5: (color online) a) Sketch of a cylindrical quantum 
dot embedded into the nanocylinder and surrounded by the 
host material. We consider that the dot yields an attractive 
potential V{r, z), represented in b) by a rectangular quantum 
well of depth Wb — —0.05eV and width w — 8nm. The radius 
of the dot is R' — Inm. 

chosen a small value for the depth of the quantum well, 
Wb = — 0.05eF, there are significant deviations in the 
tunneling coefficient from the step like characteristic, see 
Fig. joj Just before a new channel gets open, below 
there is a dip, i.e. sharp drop, in the tunneling coefficient. 
These dips are owing to modification of the tunneling co- 
efficient due to the evanescent (closed) channels^^. This 
is a multichannel effect that was until now studied only 
in Gartesian coordinates fo r quantuni w ires tailored in 
two-dimensional electron gg^^^^MIMMM^ 

The dips can be understood considering the simple 
couple mode mo del For a dot surrounded by 



host material, the scattering potential V(r, z) is not 
anymore separable, so that the scattering mixes the 
channela^^^^^^. As soon as the scattering potential is 
attractive, the diagonal coupling matrix element 



^n{r)V{r, z)(j)n{r)rdr < (39) 



acts for every channel n as an effective one-dimensional 
(ID) attractive potential^^, which always allows for at 
least one bound stat cF^^ below the threshold of the 
continuum spectrum. We have sketched in Fig. 7(a)| 
the energy spectrum of a channel n: the continuous part 
represented by continuous line is real and starts at ; 
the bound state represented by a cross (we consider for 
simplicity only one), is also real but just below the thresh- 
old. By mixing the channels, this bound state becomes a 
quasi-bound state or resonance, i.e with complex energy. 
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(a) For a channel n 



X 
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^±1 ^±2 

(b)For 2D Hamiltonian, for the first two channels 
FIG. 7: Sketch of the energy spectrum. 



whose real part gets embedded into continuum spectrum 
of the lower channel, and the imaginary part describes 
the width of the resonance. The spectrum of the 2D 
scattering problem is a superposition of the above dis- 
cussed spectra and is sketched in Fig. |7(b)| for spectra 
corresponding to channels 1 and 2. These resonances can 
be seen now as dips in the tunneling coefficient. The en- 
ergy difference between the position of the dips and the 
next subband minima E^^^ gives the quasi-bound state 
energy. The positions of the dips, i.e. the quasi-bound 
state energy, depend on the channel number n and on the 
magnetic quantum number m and, of course, on the de- 
tailed system parameters. In Cartesian coordinates the 
specific symmetry of the channels (odd and even) do not 
allow for dips in the first plateau^'^. In the cylindrical ge- 
ometry this symmetry is broken, so that we obtain a dip 
in front of every plateau. Our numerical method allows 
for high energy resolution in computing the tunneling co- 
efficient, so that we were able to find the dips also in front 
of higher-order plateaus. 

Further insight about the quasi-bound states of the 
evanescent channels can be gained looking at the wave 
functions, whose square absolute value r, z)p 

gives the localization probability density. Considering 
that the scattering states are orthonormalized in the 
general sense, the more appropriate quantity to analize 
would be the local density of states 



g{E,r) = \^i:\E,r,z)\^gsnm{E), 



(40) 



which differs from the localization probability density 
just by ID density of states. For this reason we plot 
the localization probability density in arbitrary units. 

Our numerical implementation based on the R-matrix 
formalism allows us to produce high resolution maps of 
the wave functions inside the scattering region, see Eq. 
(23). In Figs, [sjlolis represented the localization proba- 



bility density \ipn{E^ ^)P of an electron, incident from 
source (s = 1) and has a total energy corresponding to 
the dips in Fig. [6j The total energy E and the chan- 
nel n, on which the electron is incident, are specified 
at every plot. Let discuss Fig. |8(a) The total en- 
ergy E = 0.238eV is less than the energy of the sec- 
ond transversal mode, ~ 0.244eV, so that only first 
channel is open, thus the incident wave from the source 



contac t is no de-less in r-direction. But, as it can be seen 
in Fig. |8(a)[ the scattering wave function inside the scat- 
tering region has a node in the r-direction, i.e. position 
in r where the wave function is zero. This means that 
the wave function corresponds to the quasi-bound state 
splitting off from the second transversal mode, which is 
an evanescent one. The quasi-bound state is reachable 
now in a scattering formulation due to channel mixing. 
The wave function has a pronounced peak around the 
scattering potential, i.e. z G [— 4,4]nm, which decreases 
exponentially to the left and to the right. To the left 
of the scattering potential one observes the interference 
pattern produced by the incident wave and the reflected 
one, while to the right only the transmitted part exists. 

The wave function considered in Fig. |8(b)| has the 
energy less than the third transversal channel, £^^3 = 
O.GOOGeF, so that the incident part of the scattering state 
on the second mode n = 2 has one node in r-direction. 
But the scattering function shows inside the scattering 
region two nodes in the r-direction, so it corresponds to 
a quasi-bound state splitting off from the above evanes- 
cent channel, the third one. One gets similar pictures for 
all m- values, with the difference that for m 7^ the wave 
functions are zero for r = 0, like it is shown in Fig. [9] 
for the case \m I = 1. In Figs. [8] and [9] one can observe 
that the transmitted part of the scattering wave func- 
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FIG. 8: (color online) Localization probability density, 
\iljn\E,r, z)]"^ , for an electron with m = 0, incident from 
reservoir s = 1 into channel n and with total energy both 
indicated in the captions. The energies are the dips in Fig. 
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FIG. 9: (color online) Localization probability density 
|'0nm(^, r, z)!^ for an electron, with \m\ = 1, incident from 
reservoir s = 1 into channel (nm) and with total energy E, 
indicated in the captions. The energies are the dips in Fig. 
6(b)] for Wb = -OmeV. 



tion is zero, in agreement with the resonant backscatter- 
ing specific to the quasi-bound states of the evanescent 
channels 

The extension of the quasi-bound state of an evanes- 
cent channel is given outside the scattering region by 
the exponential decaying functions exp[/^in^(z + d^)] for 
z > —dz and exp[— /^2nm(^ — dz)] foi" z < dz^ where 
f^snm = ^ksnm\ and ksnm IS defined in (14). This means, 



the closer the resonance to the threshold of the evanes- 
cent channel the slower the exponential function 
decreases, yielding long exponential tails into the leads. 
This can be clea rly se en for the quasi-bound state rep- 
resented in Fig. |9(a)[ whose quasi-bound state is just 

(1) 

0.92meV below the subband minimum 

see Fig. |6(b)[ Since the localization probability density 
enters the quantum calculation of the charge distribution 
one gets difficulties in setting the correct boundaries for 
the Hartree calculations, i.e. Schrodinger-Poisson sys- 
tem. This has to be studied in a future work. 

Increasing the strength of the attractive potential to 
Wb = —0.15eV one can see more dips^^ in the tunneling 
coefficient in Fig. [6] Interesting is that there are two 
dips in the first and second plateau for m = 0, while 
for \m\ = 1 there is only one dip in every plateau. This 
can be easily understood if one thinks at the effective 
attractive potential Vnn^ Eq. (39), created for every sub- 
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FIG. 10: (color online) Localization probability density 



(1) 



E,r,z)\ for an electron with m = 0, incident from 



reservoir s = 1 into channel n and with energy E, both indi- 
cated in the captions. The energies correspond to the dips in 
Fig. [6]for Wb = -0.15eV. 
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band n. In the case of m 7^ 0, the transversal modes 
(t>n^\r) are zero on the cyhnder axis, r = 0, so that the 
effective potential for every subband is weakened. To 
confirm that the dips correspond to higher-order quasi- 



bound states, we plot in Fig. [TO]the probability density of 
the scattering states at the energies corresponding to the 
two dips in every plateau for m = . For the dips on the 
first plateau. Figs. |lQ(a)| and |lQ(b)[ both scattering wave 
functions have a node in r-direction corresponding to the 
transversal channel n = 2. For the dips on the second 
plateau. Figs. |lQ(c)| and 10(d), the wave functions have 
two nodes in r-direction corresponding to the transver- 
sal channel n = 3. But looking in the z-direction, the 
scattering state for the lower energy dip in every plateau 
is node- less, while the one for higher-energy dip in every 
plateau has a node at z = 0, which is evidence of the sec- 
ond quasi-bound states of the next evanescent channel. 



B. Core/shell quantum ring 

Now, consider the same rectangular quantum well but 
off-centered. This would corespond to a quantum ring 
embedded into the nanocylinder, as sketched in Fig. 
ll(a)| and could be realised in a core-shell heterostruc- 
ture with suplimentar structuring along the nanowire. 
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FIG. 11: (color onine) a) Sketch of a quantum ring surrounded 
by the host material. We consider that the ring yields and at- 
tractive potential V{r^z) represented in b) by an off-centered 
rectangular quantum well of depth Wh — — O.OSel/ . The thick- 
ness of the ring is Inm, between Ri — 2nm and R2 — 3nm, 
and the width of the ring is 8nm. The radius of the cylinder 
is = 5nm. 



ted in Fig. [T2j showing the characteristic dips due to 
the quasi-bound states of the evanescent channels. The 
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FIG. 12: Total tunneling coefficient as function of incident 
energy E for magnetic quantum numbers m = for a ring 
surrounded by the host material. The structure parameters 
are as in Fig. |11| The symbols show the energies, at which 
the wave functions are analyzed in the next graphs. 

localization probability density for the energies marked 



with symbols in Fig. [12] are plotted in Fig. 13 By shift 
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FIG. 13: (color online) Localization probability density 
|'0im(£^, r, z)|^ for an electron with m = 0, incident from 
reservoir s = 1 into channel n and with energy E indicated in 
the captions. The energies correspond to the symbols in Fig. 



The tunneling coefficient T'^^^E) for m = is plot- 



ing the potential from the cylinder axis and keeping the 
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same parameters as for the quantum dot surrounded by 
the host material, one can recognize the same behavior 
of the wave functions corresponding to the quasi-bound 
states as in the previous case. This means that the quasi- 
bound states of an evanescent channel extend over the 
whole width of the nanowire, independent where the scat- 
tering potential is located in the lateral direction. Similar 
results hold for m 7^ 0. Considering a deeper quantum 



Fano lin^I^I^. Both intersubband transmission probabil- 
ities T^l^ and coincide and show asymmetric Fano 
lines with zero minima. Now it is clear that there is no 
inter ference pattern to the right of the quantum ring in 
Fig. 15(a) because both transmission probabilities T^l^ 



and t|P" are zero for the second quasi-bound state. One 



recognizes in Fig. 15(a) for the first channel strong in- 




E[eV] 



FIG. 14: (color online) The tunneling coefficient T^^-* for dif- 
ferent magnetic quantum numbers m for a ring surrounded 
by the host material. The parameters are as in Fig. [11] but 
Wb — — 0.15ey. The symbols show the energies, at which the 
wave functions are analyzed in the next graphs. 

well one is surprised to see in Fig. [14] that for m = only 
one dip appears in the first plateau but two dips in the 
second plateau. This can be understood considering that 




{h)E = 0.599ey, n = 2 



FIG. 15: (color online) Localization probability density 



• r .1 . 11 1 \ u J ^ rr \^ri (E.r, z)\ for an electron with m = 0, incident from left 

II the transversal channels 0n (r) has a node at the on- y ) ]^ 1 1 / x 1 o /i x 1 ..1 



centered position of the scattering potential, then Vnnj 
Eq. (|39| , is being weakened allowing for less quasi-bound 
states. Interesting is to take a closer look at the scatter- 
ing states corresponding to the second quasi-bound state 



into both open channels n = 1 (a) and n = 2 (b) and with 
energy E corresponding to the second quasi-bound state on 
the second plateau in Fig. |14| 



on the second plateau, marked with a symbol in Fig. 14 
For this energy, there are two open channels, and we have 
represented in Fig. [15] the localization probability den- 
sity for both of them. One can recognize immediately 
the structure of the wave function with two nodes in r- 
direction corresponding to the third evanescent channel 
and one node in z-direction specific to the second quasi- 
bound state. Specific for the quasi-bound states of an 
evanescent channel is also the exponential decaying far 
from the scattering potential. But the interference pat- 
terns on the left and on the right of the scattering po- 
tential are quite different for these two scattering wave 
functions. This can be explained by looking at intrasub- 
band and intersubband transmission probabilities repre- 
sented in Fig. [16] and which give detailed information 
about channel mixing. The intrasubband transmission 
T^j^ is stronger influenced by channel mixing, showing 
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FIG. 16: (color online) Intrasubband and intersubband trans- 
mission probabilities for a quantum ring (Fig. 



Wt 



-0.15eV and for m = 0. The curves for T. 



(1) 



11) with 



and T^l^ 



two pronounced dips, while the intrasubband transmis- coincide. The vertical doted line shows the third subband 
sion T22^ shows only one dip and a second asymmetric minimum ^^3. 
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terference pattern between the incident part and the re- 
flected part. For the scattering wave function incident on 
the second channel T22^ has values close to 1, and also 



T^j^ h as a maximum. In such a way one sees in Fig. 



15(b) right to the scatterer an interference pattern be- 
tween the transmitted wave in the first channel and the 
one transmitted in the second channel. One can recog- 
nize far from the scattering potential the structure of the 
second channel with a node in the r-direction. 



C. 



Double-barrier heterostructure along the 
nanowire 



In Fig. |17(a)| is sketched a double-barrier heterostruc- 
ture along the cylindrical nanowire. Such systems with 
sharp interfaces between the layers are realized exper- 
imentally based on InAs/Inl^ or on GaAs/AlGaAP. 
We consider rectangular barriers of height V5 = O.beV^ 
widths b = 4nm, and the width of the rectangular quan- 
tum well w = 8nm, as is plotted in Fig. 17(b) The 
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FIG. 17: (color online) a) The sketch of a double-barrier het- 
erostructure along the nanowire and b) the scattering poten- 
tial V{r,z). The height of the barriers is Vb = 0.5eV, the 
widths of the barriers is 6 = 4nm, the width of the quantum 
well is w = 8nm and the radius of the nanowire is = 5nm. 

total tunneling coefficient T'^^^E) for m = is plotted 
in Fig. 18(a) in linear scale. The barriers supprese the 
transmission, except for a series of sharp peaks due to 
the quasi-bound states between the barriers. With verti- 
cal dashed lines we have represented the energies of the 

,(0) 



can observe that the tunneling coefficient can reach val- 
ues higher than 1, if there are more channels open. In 
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(b) logarithmic scale 

FIG. 18: (color online) Total tunneling coefficient in lin- 
ear scale (a) and logarithmic scale (b), as function of inci- 
dent energy E for a double-barrier heterostructure along the 
nanowire, as depicted in Fig. |17| The symbols represents 
the poles: their real part on x-axis and the imaginary part on 
right y-axis. The peaks are indexed by (n, i), where n denotes 
the incident channel and i denotes the resonance between the 
barriers. The same indexes are used in Tabled 

case of no applied bias between source and drain con- 
tact, the scattering potential is separable and has vari- 
ations only in z-direction V{r^z) = V{z), where V{z) 
describes an ID double-barrier potential. The scatter- 
ing does not mix the channels, so that the total tun- 
neling coefficient is given by summation of the intrasub- 
band transmission probabilities for every open channel 
T^^\E) = ^j^Tnn{E). The intrasubband transmission 
probability on every open channel is the transmission 
through a double-barrier structure, but shifted with the 
transversal energy of the channel E^^^ . So that it can be 



ffist three transversal channels £^x,n' ^ = 1,2,3. One also computed as for a double-barrier resonant tunneling 
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TABLE I: Localization probability density \ipn\Eres,i,r, z)\^ 
for an electron with m = incident from reservoir s — 1 into 
channel n and corresponding to the resonance i between the 
barriers. The axis of abscissae is z G [—8, 8]nm, and the axis 
of ordinates is r G [0, 5]nm for all plots. 



(DBRT) diode, t£\e) = TDBRT{e + eJ^)- 

This identity can be used as a verification of the nu- 
merical implementation of our method, because the first 
quantity is computed with the 2D code, while the second 
quantity is computed with the ID cod^l. This is ex- 
plicitly illustrated in a lo garithmic plot of the tunneling 
coefficient in Fig. |18(b)[ We represent here by vertical 
dashed lines the positions of the first three transversal 
channels ^^n' ^ ~ 1,2,3. By dot-dashed lines we have 
represented the tunneling coefficient through the DBRT, 
but the energies are shifted with the transversal chan- 
nel energy, as is written in the legend. The curve for 
the first channel (red dotted line) is just under the to- 
tal transmission curve (black continuous line). One can 
observe that the ID double-barrier potential allows for 
four resonances (quasi-bound states) between the barri- 
ers, i.e. every dashed curve has four peaks below the 
barrier height. 

On the same plot we have plotted with symbols the 
poles of the current scattering matrix S, computed using 
the method presented in Ref. [13] recently developed for 
2D geometries^^ . The real part is on the axis of abscissae, 
while the imaginary part is on the right axis of ordinates. 
One can see that among the poles there exist resonant 
ones, marked by filled symbols, with very low widths, i.e. 
r < 10~^eV^ and which are very well separated from the 
others. Using the same axis of abscissae for the real part 
and for the tunneling coefficient, one can directly see that 
to every resonant pole corresponds a transmission peak. 
Increasing the energy, but keeping the same channel, the 
widths of the poles increase, and so the widths of the 
transmission peaks. 

This physical interpretation of the tunneling coefficient 



peaks allows us to label them in Fig. 18 by a pair of num- 
bers (n,i), where n describes the incident channel and i 
describes the resonance (quasi-bound state) between the 
barriers. 

In Table [T] are represented the localization probabil- 
ity densities \ipn\Eres,i^ for the double-barriers re- 



gion, z e [— 8,8]nm and r G [0,5]nm, for every indexed 
peak in Fig. [18] One can observe that the wave func- 
tions have pronounced maxima between the barriers and 
decrease very quickly inside the barriers. They are local- 
ized between the barriers, corresponding, indeed, to reso- 
nances (quasi-bound states) between barriers and not to 
quasi-bound states of evanescent channels. Furthermore, 
the order i of the resonance between the barriers gives 
the number of nodes in the z-direction, namely i — 1, 
while the channel number n gives the number of nodes 
in the r-direction, namely n — 1. This provides a picture 
of the orbitals of the "artificial atom", which represents 
this quantum structure^^. 
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FIG. 19: (color online) Transmission coefficient T^^-* for |m| = 
1 vs. total energy E for a double-barrier heterostructure along 
the nanowire as depicted in Fig. |17| The peaks are indexed 
by (n, i), where n denotes the incident channel and i denotes 
the resonance between the barriers. The same indexes are 
used in Table [n] 
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TABLE II: (color online) Localization probability density 
for an electron with |m| = 1 incident from 
reservoir s = 1 into channel n and corresponding to the 
resonance i between the barriers. The axis of abscissae is 
z G [—8, 8]nm, and the axis of ordinates is r G [0, 5]nm for all 
plots. 
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Similar behavior is observed for higher magnetic quan- 



tum numbers m. In Fig. [19] is represented the total tun- 
neling coefficient and in Table [TI| the localization proba- 
bility densities for the indexed peaks for the case |m| = 1. 
The positions of the transmission peaks vary for different 
m- values due to the dependence on m of the transversal 
energy channels The scattering wave functions at 

resonances for different m-values have different positions 
of nodes in the r-direction, and for any m 7^ they are 
zero on the cylinder axis. 



IV. SUMMARY AND DISCUSSION 

We have presented a general theory for computing the 
scattering matrix and the scattering wave functions for a 
general finite-range extended scattering potential inside 
a cylindrical nanowire. This formalism was applied to a 
variety of model systems, like a quantum dot, a quan- 
tum ring and a double-barrier heterostructure embedded 
into the nanocylinder. We have recovered the features 
for a nonseparable attractive scattering potential in a 
multi-channel two-probe nanowire taylored in the two- 
dimensional electron gas. The difference to the Carte- 
sian geometry is that every magnetic quantum number 
defines a two-dimensional scattering problem with differ- 
ent structure of dips for the same scattering potential. 
How many of these problems have to be solved depends 
on the further physical quantity calculation. Further- 



more, the cylindrical symmetry does not yield the same 
"selection rules" for tunneling coefficient as the Carte- 
sian symmetry, so that dips could be observed in every 
subband. For stronger attractive potential more than 
one dip can appear due to the higher-order quasi-bound 
states of the next evanescent channel. For quasi-bound 
states localized between barriers, it was possible to com- 
pute the poles of the scattering matrix, which provide a 
quantitative characterization of the resonances. Further- 
more, the peaks of resonant tunneling can be indexed by 
channel number and resonance index. Detailed maps of 
localization probability density sustain the physical in- 
terpretation of the resonances (dips and peaks) found in 
the studied nanowire heterostructures. 

It will be the subject of next works to see, how the 
buildup of charge around the scattering nonseparable at- 
tractive potential influences the overall electrical charac- 
teristics of the nanowire-based devices. 
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